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We present the threshold resummation for the single spin asymmetry in heavy flavor production 
by linearly polarized photons. We analyze the soft-gluon contributions to the asymmetry at flxed 
order in Us to the next-to-leading logarithmic accuracy. Our analysis shows that, contrary to the 
production cross section, the azimuthal asymmetry is practically insensitive to soft radiation. Our 
calculations of the asymmetry up to the 6th order in Us lead only to small corrections (of order of 
few percent) to the Born predictions at energies of the fixed target experiments. Fast convergence 
of the perturbative series for the azimuthal asymmetry is due to the factorization properties of 
the photon-hadron cross section. We conclude that measurements of the single spin asymmetry 
would provide an excellent test of pQCD applicability to heavy flavor production. 
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(D ■ 

CSJ ; I. INTRODUCTION 

, Presently, the basic spin-averaged characteristics of heavy flavor hadro-, photo- and electroproduction are known 
^ ' exactly up to the next-to- leading order (NLO) (see for a review). Two main results of the explicit pQCD 
. calculations can be formulated as follows. First, the NLO corrections are large; they increase the leading order 
(LO) predictions for both charm and bottom production cross sections approximately by a factor of 2. For this 

• reason, one could expect that the higher order corrections as well as the nonperturbative contributions can be 
, essential in these processes, especially for the c-quark case. Second, the one-loop predictions are very sensitive to 

• standard uncertainties in the input QCD parameters. In fact, the total uncertainties associated with the unknown 
values of the heavy quark mass, m, the factorization and renormalization scales, fip and /i^, A^qcd and the parton 
distribution functions are so large that one can only estimate the order of magnitude of the NLO predictions for 

^pH. total cross sections [||J^. For this reason, it is very difficult to compare directly, without additional assumptions, the 
O il NLO predictions for spin-averaged cross sections with experimental data and thereby to test the pQCD applicability 
^ • to heavy quark production. 

During the recent years, the role of higher order corrections has been extensively investigated in the framework 
of the soft gluon resummation formalism. For a review see Ref. Q. Soft gluon (or threshold) resummation is 
. based on the factorization properties of the cross section near the partonic threshold and makes it possible to 
5— I ' resum to all orders in as the leading (Sudakov double) logarithms (LL) and the next-to- leading ones (NLL) 

Formally resummed cross sections are ill-defined due to the Landau pole contribution, and some prescription must 
be implemented to avoid the renormalon ambiguities. To obtain numerical predictions for the physical cross sections, 
a few prescriptions have been proposed [^-p^. Physically, the choice of a resummation prescription implies the 
introduction of an effective scale, separating the hard (perturbative) and soft (nonperturbative) gluon contributions. 
Another open question, also closely related to the convergence of the perturbative series, is the role of the subleading 
logarithms which are not, in principle, under control of the resummation procedure [ pT| , p^ . Numerically, the higher 
order corrections to the heavy quark production cross sections can depend significantly on the choice of resummation 
prescription [p^ . 

Since the spin-averaged characteristics of heavy flavor production are not well defined quantitatively in pQCD it 
is of special interest to study those spin-dependent observables which are stable under variations of input parameters 
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of the theory |13[ . In this paper we analyze the charm and bottom production by hnearly polarized photons, namely 
the reactions 

-f^+h^Q + X[Q]. (1.1) 

We consider the single spin asymmetry parameter, A(S), which measures the parallel-perpendicular asymmetry 
in the quark azimuthal distribution: 

^iS,^) = ^2|Ai(i + A(5)7,^cos2^), (1.2) 

where 

Here cr^)J^(S') is the unpolarized cross section, da{S, ip) = ^{S, ip) , Vj is the degree of linear polarization of the 

incident photon beam, \/S is the centre of mass energy of the process (1.1) and (p is the angle between the beam 
polarization direction and the observed quark transverse momentum. 



The Born level predictions for the single spin asymmetry in (1.1) as well as the contributions of nonperturbative 
effects (such as the gluon transverse motion in the target and the heavy quark fragmentation) have been considered 
m H]. The following remarkable properties of A{S) have been observed: 



The azimuthal asymmetry (1.3) is of leading twist; in a wide region of initial energy, it is predicted to be 



about 20% for both charm and bottom quark production. 

• At energies sufficiently above the production threshold, the LO predictions for A(S) are insensitive (to within 
few percent) to uncertainties in the QCD input parameters: m, ^n, ^p, Aqcd Si,iid in the gluon distribution 
function. This implies that theoretical uncertainties in the spin-dependent and spin-averaged cross sections 



(the numerator and denominator of the fraction (1.3), respectively) cancel each other with a good accuracy. 



• Nonperturbative corrections to the 6-quark azimuthal asymmetry A{S) due to the gluon transverse motion in 
the target are negligible. Because of the smallness of the c-quark mass, the analogous corrections to A{s) in 
the charm case are larger; they are of the order of 20%. 



In the present paper, the soft-gluon corrections to the asymmetry ( |1.3| ) are investigated. Using the methods of 
Rcfs. 1^-0] for the resummation of soft gluons, we resum the Sudakov logarithms for the polarized cross section 



(1.2) in moment space to the NLL accuracy. Considering the obtained expression as a generating functional of 
the perturbative theory, we re-expand it in as in momentum space and analyze the fixed order predictions for the 
azimuthal asymmetry. Our main results can be formulated as follows: 



Contrary to the production cross sections, the asymmerty (1.3) in azimuthal distributions of both charm and 



bottom quark is practically insensitive to radiative corrections at fixed target energies. This implies that large 



soft-gluon contributions to the spin-dependent and unpolarized cross sections cancel each other in ( L3 ) with 
a good accuracy. 

At the NLL level, the NLO and NNLO predictions for A{S) affect the LO resuhs by less than 1% and 2%, 
respectively. 

Our computations of higher order corrections indicate a fast convergence of the perturbative series for A{S). 
We have calculated the asymmetry parameter A{S) up to the 6th order (N^LO) in as to the NNL accuracy 
and have found that corresponding corrections to the Born level result are of order of few percent. 

Our analysis shows that high perturbative stability of A[S) is due to the factorization properties of the photon- 
hadron cross section. We conclude that, in contrast with the production cross sections, the single spin asymmetry in 
heavy flavor photoproduction is an observable quantitatively well defined in pQCD: it is stable, both parametricaly 
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and perturbatively, and insensitive to nonperturbative corrections. Measurements of the azimuthal asymmetry in 
bottom photoproduction would provide an ideal test of the conventional parton model based on pQCD. 

Concerning the experimental aspects, the azimuthal asymmetry in charm photoproduction can be measured at 
SLAC where a coherent bremsstrahlung beam of linearly polarized photons with energies up to 40 GeV will be 
available soon |l^Jl5[] . Due to the c-quark low mass, data on the _D- meson azimuthal distributions would make it 
possible to clarify the role of subleading twist contributions |13|. 

The paper is organized as follows. Section II contains the derivation of the rcsummcd formula for the spin- 
dependent photoproduction of heavy flavor in the single-particle inclusive kinematics. In Section III we analyze the 
NLO and NNLO predictions for the azimuthal asymmetry. We check the quality of the NNL approximation against 
available explicit results and discuss the subleading logarithms contribution. Higher order corrections and the role 
of the gluon distribution function in perturbative stability of the asymmetry are considered in Section IV. 



II. SOFT GLUON RESUMMATION FOR POLARIZED CROSS SECTION 

In this Section we carry out the resummation of the Sudakov logarithms for the spin-dependent cross section of the 
reaction 

7(^) T +Hkh) ^ QIpq) + xmpx) (2.1) 

to next-to-leading logarithmic (NLL) accuracy to all orders of the perturbative expansion. In the single-particle 
inclusive (IPI) kinematics, the overall invariants are defined as 

S ^ {k^ + kh f ; Ti = (kfi - pq f - m^; 

Si = S + Ti + Ui; Ui = [k^ ~ pq f - m\ (2.2) 

At the parton level, the dominant subprocess is the photon-gluon fusion: 

7(fc^) + gikg) QIpq) + X[Q]ipx), (2.3) 
where the corresponding kinematical variables are 

s = {k^ + kgf' = zS; ti = {kg - pqf' - = zTi; 

S4 = s + ti+ui; ui^Ui, (2.4) 

with kg — zkh, m the heavy quark mass, while 54 measures the inelasticity of the partonic reaction. At the Born 
level, 0{aemO:s), the only partonic subprocess which is responsible for heavy quark photoproduction is the two-body 
photon-gluon fusion: 

l{k^)+9{kg)^Q{pQ) + Q{p^). (2.5) 

To take into account the NLO contributions, one needs to calculate the virtual 0{aemC(f) corrections to the Born 
process ( pj| ) and the real gluon emission: 

7(^7) + 9{kg) ^ QiPQ) + QiPg) + 9{Pa)- (2.6) 

We neglect the photon-(anti)quark fusion subprocesses as well as the so-called hadronic or resolved component 
of the photon. This is justified as their contributions vanish at LO and are small at NLO at the energies under 
consideration 

The factorized single-heavy quark inclusive cross section for photon- hadron collisions, Aa^h, has the form of a 
convolution of the perturbative short-distance cross section, Aa^g, with the universal parton distribution function 
't>g/h- 
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da^h {S, Ti, Ui,ip, fiF,as{fJ.R)) = / dz(f)g/h{z, HF)da^g (s, ii, ui, Lp, ^if, as(Mlj)) 



(2.7) 

where hf and fin are the factorization and renormaUzation scales, respectively. Note that here and in the following 
da^fi and da^g denote any relevant (^-dependent differential distribution. Replacing the incoming hadron by the 
gluon and taking the Laplace moments, the above convolution simplifies to a product: 



dajg{N,ip) = (j)g/g{Nu) da^g{N,ip). 



The moments for dcr-^g are defined by 



da^g{N, if) 



dsA 



.e-^s,/m-da^gisi,v), 



(2.8) 



(2.9) 



with N the moment variable. The upper limit of this integral is not important for large N and may be put at 1. 
Similarly, the moments for 4>g/g with respect to z have the form 



4>g/g{Nu) = 



dze" 



-Af„z 



4>g/g{z), Nu = N{-ui/m'^), 



where definition of is given for the IPI kinematics p3||lq|. 

The short-distance perturbative cross section da-yg still be sensitive to the collinear gluon emission, pg^x 
separate these collinear effects from the hard scattering, a refactorization is introduced p7|,p[: 



dajg{N,ip) = tpg/g{Nu)Hjg{ip)Sjg 



( m 



(2.10) 



0. To 



(2.11) 



where tpg/g is the center-of-mass parton distribution that absorb the universal collinear singularities associated 
with the initial-state gluon while S^g is the soft-gluon function that describes the non-coUinear soft gluon emission. 
The mass of the heavy quarks protects the final state from collinear singularities. The hard-scattering part of 
cross section, H^g, is free of soft-gluon effects and thus independent of TV. In the photoproduction case, both H-yg 
and S-^fg are simply functions, and not matrices in color space, in contrast with heavy quark or jet production in 
hadron- hadron collisions [Q. 



Comparing Eqs. (2.8) and (2.11), we derive 



da^g{N, ip) 



V'g/s(^")a f m \ 

~ ij' 



[NfiFj 



H^gip). 



(2.12) 



The functions i^g/gl4>g/g and S-yg originate from the collinear, pgj^ — > 0, and soft, Pg — > 0, limits. Since the 
azimuthal angle tp is the same for both 717 and QQ center-of-mass systems in these limits, only the hard function. 



Hyg, is (f- dependent in the right-hand part of (2.12). 

In the MS factorization scheme, at NLO and to NLL accuracy, the ratio "^g/g/^g/g is [p^ 



-~ [J^u, fJ.F) = H 



\n^Nu+{l + In 



(2.13) 



where = NuQ*'^ with 7^ the Euler constant. The function il>g/gl4'g/g summarizes all leading logarithms of N 
and a part of the next-to- leading ones. (At NLO, they arc In^ iV and IniV, respectively). 5-yg-function contributes 
only at the next-to-leading level of In N . 

To resum the leading, O (^a^ In^" , and next-to-leading, O {a'^lv?"^^^ N^, logarithms for da^g (iV, </?) to all 

orders of n, we use the methods of Refs. (^-0] based on the solving of the appropriate evolution equations for each 
of the functions in (^.12 ). The formal result for the resummed spin-dependent cross section in moment space is: 
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/■™/^ An' 

/ ^2Rer-K(M'^)) 

Jm r 



exp 



(2.14) 



The first exponent in (2.14) resums the iV-dependence of the ratio '4'g/g/<t>gig and is given in the MS scheme by 



E,{Nu) = I ^ (1 



g/a/'Pa/a 
1 

— {as (/i'^)) + -i^g {as (w^m^)) 

J 2 ^ 



(2.15) 



The function Ag is known at two loops 0JT|, (a^) = a^CA/Tr + (a^/Tr)^ CaK/2, K ^ Ca (67/18 - ttVO) - 
5nf/9, where is the number of active quark flavors, and Vg = ^UsCaIt^ 0- 

The second exponent controls the evolution of the soft function from scale m/N to m and is given in terms of 
r^^, the soft anomalous dimension. F^^ was calculated at one loop in JTot : 



n'{o.s) 



2ti 



In- 



u^t^ 



l-|f)(l + .,a) 



Here 



L, 



1 - 2m2/s 



In 



1 + /3 



^ = v/l-4m2/s, 



(2.16) 



(2.17) 



while Ca = iVc and Ci? = (7V^ — l)/(2A^c), where Nc is the number of colors. The product H^gS^g on the first line 



of Eq.(2.14) is determined from matching at lowest order in as to the Born result. 



The exponent in Eqs.(2.14), (2.15) can only be interpreted in a formal sense since the corresponding integration 
over Lo near Aqcd I'm is not defined unambiguously due to the Landau pole in the coupling strength as ■ To avoid 
this soft gluon divergence, several prescriptions have been proposed |p[-JlO|. However, if one expands the exponents 
in the resummed cross section at fixed order in a^, no divergences associated with the Landau pole are encountered. 



III. NLO AND NNLO PREDICTIONS 
A. Partonic Cross Sections 



In this section we expand the resummed cross section (2.14) to one and two loop order and invert back to momentum 
space using the equation |l^ 



(- In iv)'^^ = (; + 1) / dsie-^'^' 



Si 



O (In^^iV 



where IPI singular "plus" distributions are defined by 



W (s4/m2 

S4 



lim 



In' (s4/m2) 



1 



(^4-e) + r^ln'+M32)'5(54) 



Si 



l + l 



For any sufficiently regular test function /i(s4), Eq.(3.2) gives 





{si/m'^y 




/ (lsih[si) 


.-I 


Jo 


S4 





(3.1) 



(3.2) 



ds4 [h{s^) - h{0)] ^ J_^Q) i^i+i (s--/to2) . (3.3) 



S4 



Eq.( |3.lD provides us the NLL approximation for both spin-dependent, dAtr^g, and unpolarized, da^g, differential 
cross sections: 



5 



S 111 I 11 11 ' T I , 

(\ipQ.t\Au\ 27r \ d^idui diiditi 



(3.4) 



where V.^ is the degree of the photon beam polarization; kp is the angle between the observed quark transverse 
momentum, pq^t, and the beam polarization direction. 

To NLL accuracy, the perturbative expansion for the partonic cross section can be written in a factorizcd form 

as 



2d2(A)(T^g (s,ii,ui) 



^^^^^^ (A)B^;"(s,ti,Ki)|<5(5 + ti+tii) + ^^(^^j 

with the Born level hard parts {A)B^°"' given by ||,|3| 

ti ui Arn^s 



(3.5) 



(s,ii,ui) = 27reQaemasTf 



Ul tl 



t-iU 



I'Ul 



1 - 



(3.6) 



Ai3^°™ (s, ii, Ul) = 27reQQ;emas?F 



4m s 



1 



(3.7) 



In (|j) and (||3), Tp =Tr{T''T'')/{N^ - 1) 
electromagnetic coupling constant. 

The NLO and NNLO soft gluon corrections to NLL accuracy are 



tiUi \ tiUi ^ 

1/2 is the color factor, and eg is the quark charge in units of 



In [si/m? 

Si 



S (S4) In 





■ 1 " 























ln(^ 



2Cf 

c7 



(l + ReL;3)+ln ^ 



and 



In^ (s4/m2) 

S4 

In^ (s4/m^) 

S4 

In [si/m?' 



S4 



In ( ^ 1 - 



1 + ln 



1 +ln 



- 1- 



2Cf 

c7 



-vri-^)(i 



(l + ReL,,) 
(1 + KcLfi) -[ 



2b2_ 

3Ca 



In 



In 



-Ul 



Ca 



-In 
2 



ln^(^ 



In 



1_&2_ 

2Ca 



where we use /i — — M-Rj and 62 is the first coefficient of the [3 (Q;s)-function expansion: 

dlnofs 



/9(«s) 



d In /i^ 



^6fe+l(a./7^)^ 



(3.8) 



(3.9) 



(3.10) 



62 = (UCa - 2n/) /12, 63 = (34Ci - IOCa^/ - leC^n/) /48. In (gj) and (|3J), we have preserved the NLL 
terms for all the orders of the scale-dependent logarithm In*^ (/z^/m^), (A: = 0, 1, 2). We have checked that the result 
( |3.8[) agrees to NLL accuracy with the exact 0{aemOi^) calculations of the unpolarized photon-gluon fusion given 
in |22|p^ , and that the Eqs.(3^) and ( |3.9| ) coincide completely with the corresponding ones obtained in (l||l| for 
the clcctroproduction case. 
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To perform a numerical investigation of the results (3.8) and (3.9), it is convenient to introduce for the fully 
inclusive (integrated over ti and ui) cross sections, a-^g and Aa^g, 



/•(l+/3)s/2 

(A)a^,(s) = / d(-ii) 



max 
S4 



S + ti + 



ds4 



d^(A)a,g 



(s,ti, 54), 



the dimensionless coefficient functions (A)c'-'^''^ 



(A)0-7g(?7,M ) 



^(4™,(^^))^^(A)c('=^')(^)lh 



fc=0 



(=0 



1, 



(3.11) 
(3.12) 

(3.13) 



where the variable rj measures the distance to the partonic threshold. 

Concerning the scale- independent coefficient functions, only c^^'°-' is known exactly [^. As to the /^-dependent 
coefficients, some of them can by calculated explicitly using the renormalization group equation: 



d(A)a-;(s,Ai) 



dz(A)a:-(zs,M)Pg<,(^), 



(3.14) 



d In /i^ J p 

where p = Am? j s, {/\)a''^g{s, p,) are the cross sections resummed to all orders in as and Pgg{z) is the corresponding 
(resummed) Altarelli-Parisi gluon-gluon splitting function. Expanding Eq.(3.14) in a^, we find 

1 



(A)c' 



(1,1) 



47r2 



62(A)c("'°)-(A)c(°-")0Pg(o) 



(A), 



.(2,1) _ 



(47r2 



63(A)c("'0)-(A)c(°'°)®P« 



47r2 



262(A)c(i-0)-(A)c(i^°)®pW 



(3.15) 
(3.16) 



(A), 



.(2,2) _ 



47r2 

1 



(47r2)' 

where the convolutions (A)c*^'^''-' (X 



6.(A)c(i^i)-l(A)c(i^i)«p(°)' 

62(A)c("'°) - ^62(A)c("^") ® PW + i(A)c(«^o) P^ ^ 



Pgl' are defined as 



(A)c('=-')®P(^,)1(.)^ / dz(A)c('=-')(zs)P«)(z) 



with the one- and two-loop gluon sphtting functions, Pgg{z) — {as/'r:)Pg1\z) + {as/T:YPgg' {z) + given in |2 



(3.17) 



(3.18) 



Note that Eq.( 3.15| ) agrees with the exact result for c'^^'^-' given in |2j 



With Eqs.(|3.15D-(p7|) in hand, we are able to check the quality of the NLL approximation at both NLO and 
NNLO against exact answers. In Figs.l and 2 we plot the functions €^'^'''^(77) and Ac'*^''' (77), respectively. Predictions 
of the NLL approximation (3_^), (3^) are given by dashed curves. The available exact results are given by solid 
lines. One can see a reasonable agreement up to energies n « 2. 

0, 1, 2. The NNLO curve is given only up to 



In Fig. 3 we plot the NLL predictions for the ratios ^f„^o} iv): 



77 = 2 since, at larger energies, the coefficient function c^2,o) changes sign and the ratio strongly oscillates. 

Fig. 3 shows sizable deviations of the NLO and NNLO results from the Born level ones. This is due to the fact that 
the physical soft-gluon corrections ( ^.ll| ) are determined by a convolution of the Born cross section with the Sudakov 
logarithms which, apart from factorized (5(s4)-terms, contain also non-factorable ones (see Eq.(3.3)). Kinematically, 
large values of 77 allow Si/m? > 1 that leads to significant non- factorable corrections. In other words, the coUinear 
brcmsstrahlung carries away a large part of the initial energy. Since the spin-dependent and unpolarized Born level 
partonic cross sections have different energy behavior, the soft radiation has different impact on these quantities. 
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FIG. 3. NLL predictions for the ratio („ o) iv)- Plotted are the LO (solid curve), NLO (dashed curve) and NNLO (dotted 
curve) results. 



B. Hadron Level Results 



Let us now analyze the impact of the approximate NLO and NNLO perturbative corrections on the azimuthal 
asymmetry, A{S), at hadron level. Unless otherwise stated, the CTEQ5M ||2^ parametrization of the gluon dis- 
tribution function is used. The default values of the charm and bottom mass are nic = 1.5 GeV and mf, = 4.75 



GeV. For our analysis at NLO (NNLO) we use the two-loop (three- loop 



expression for as] A4 = 300 MeV 



and A5 = 200 MeV. The default values of the factorization scale iip chosen for the A{S) asymmetry calculation 
are fip | Charm = 2mc for the case of charm production and | Bottom = "T-b ioi the bottom case Jl|,p9t. For the 
renormalization scale, /x_r, we use hr = fip- 

Our results for the single spin asymmetry A{S) in charm and bottom photoproduction at fixed target energies 
are presented in Fig. 4. For comparison, we plot in Fig. 5 the so-called X-factors for unpolarized cross sections: 



a. 



{S)/ aLO(5) and K^^^iS) = a 



_ ^NNLO 



7/1 



NLO 



(S). One can see from Figs. 4 and 5 that large 



soft-gluon corrections to the production cross sections practically (to within 1-2 percent) do not affect the Born 
predictions for A{S) at both NLO and NNLO. 
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FIG. 4. Single spin asymmetry, AiE^,), in 6- and c-quark production at LO (solid curve), NLO (dashed curve) and NNLO 
(dotted curve) to NLL accuracy. 
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FIG. 5. Energy behavior of the K{E^)-ia.ctois, k!:^)!{S) = (7^h°{S)/ cr:^°(5) (dashed curve) and k':^^{S) = a^^^°{S)/ 
a^^'^{S) (dotted curve), in b- and c-quark production at NLL level. 



Our calculations given in Figs. 3-5 represent the central result of this paper. At first sight, the situation seems 
paradoxial: large soft-gluon corrections display a strong (^-dependence in the case of the photon-gluon fusion and, 
simultaneously, are practically i^j-independent at the hadron level. A qualitative explanation of this fact will be 
given in Section IV. We shall see that sufficiently soft gluon distribution function leads to a factorization of the 
photon-hadron cross section and, as a consequence, to fast convergence of the perturbative series for A{S). 

Another remarkable property of the azimuthal asymmetry closely related to fast perturbative convergence is its 
parametric stability]^. As it was shown in |p^ , the LO predictions for A{S) are insensitive (to within few percent) 
to standard theoretical uncertainties in the QCD input parameters: m, /ij?, fip, Aqcd and in the gluon distribution 
function. We have verified that the same situation takes place at higher orders too. In particular, all the CTEQ5 
versions of the gluon density as well as the MRST parametrizations |27| lead to asymmetry predictions which 
coincide with each other with accuracy better than 1.5%. 

Note also the scaling behavior of the azimuthal asymmetry: with a good accuracy the quantity A{S) is a function 
of the only variable ph, ph = Am? /S, so that |l^ 

A[S) |Bottom~ A[Sml/ml) |charm ■ (3.19) 

This property of the asymmetry reflects its independence from Aqcd- 

Let us now discuss the hard (pg,r 7^ 0) and virtual NLO contributions, i.e. radiative corrections to the hard 
functions, H^g{s) and AHjg{s). Since these contributions do not contain the Sudakov logarithms, they are expected 
to be small near the threshold. One can see from Fig.l that it is really the case. Soft radiation describes very well 
the exact NLO results on the unpolarized photon-gluon fusion at partonic energies up to 77 « 2. Since the gluon 
distribution function supports just the threshold region, the soft-gluon contribution dominates the photon-hadron 
cross section approximately up to S/irr? ~ 10. Using the exact expression for the 717 cross section p^ , we have 
verified that the contribution originating from the region 77 > 2 makes only 1-2% from the NLO predictions for 
the unpolarized bottom production in 7p collisions at E-y < 1 TeV. In other words, radiative corrections to the 
unpolarized hard function affect the NLL predictions for the asymmetry A{S) by about -1-2% at S/Arr? ~ 10 and 
by less than 1% at lower energies. 

Presently, the exact NLO calculations of the (^-dependent cross section of heavy flavor production are not com- 
pleted |pl|. However we can be sure that, at energies not so far from the production threshold, the soft radiation 



^Of course, parametric stability of the fixed order results does not imply a fast convergence of the corresponding series. 
However, a fast convergent series must be parametricaly stable. In particular, it must be pn- and /iF-independent. 



11 



is the dominant perturbative mechanism in the polarized case too. First, the LO predictions for the (/^-dependent 
cross section are large and the Sudakov logarithms have universal, spin-independent structure. For this reason, the 
polarized heavy quark production has also a strong threshold enhancement. Second, our analysis of the exact scale- 
dependent polarized cross sections given in Fig. 2 confirms with a good accuracy the dominance of the soft-gluon 
contribution. Third, we have verified that soft radiation is a good approximation to the threshold behavior of the 
explicit NLO results on heavy quark production in 717 collisions with longitudinally polarized initial states too. 
These facts argue that radiative corrections to the (^-dependent hard function can not also affect significantly the 
soft-gluon predictions for the asymmetry at the energies under consideration. 

At very high energies, S/Am? > 10^, the dominant NLO production mechanism is the so-called flavor excitation 
(FE). This mechanism arises from the diagrams with the t-channcl gluon exchange and leads to a constant value for 
the unpolarized photon-gluon fusion cross section as s — > 00 ]2^ , p^ . In this limit, the hard component associated 
with the FE mechanism is large and can, in principle, affect the LO predictions for the asymmetry. One can expect 
that hard corrections will result in some dilution of the azimuthal asymmetry at superhigh energies because of 
different asymptotical behavior of the spin-dependent and unpolarized cross sections: at s —> 00, Cp-^\s) cx const 
while Acpg''^(s) (x Ins/s. In particular, hard corrections can be sizable in the case of charm production at ^ 1 
TeV. For this reason, we consider only the energy region up to S/Am? ~ 10. 

In the conclusion of this section, a remark about the next-to-next-to-leading logarithms (NNLL) contribution. In 
the unpolarized case, it is possible to estimate the contribution of NNLL at NNLO using the method proposed in 



Refs. [B(3,n2|. For this purpose, we can add to the NLL approximation (3.8) the subleading 5(s4)-term containing 



the exact virtual plus soft corrections Exponentiating the obtained expression, one can determine the 

coefficients of all the powers of the threshold logarithms [In''" (§4/777,^) /•S4]+ at NNLO. (Remember that the NLL 
expansion (3^) does not give all the terms with fc = 0, 1 but only those ones involving the scale ^). Since the 
exact expression for Ac*-^'"' is not presently available, in the polarized case we can analyze numerically only the 
factorized contribution of the so-called Coulomb singularity which originates from the final state interaction between 
the massive quarks and dominates in the region very close to the threshold ||2|] . We have verified that the Coulomb 
corrections to both Ac'^^'^-' and c^^'^-' are of the order of few percent and that their contribution to A{S) is negligible. 

Beyond the NNLO, not all of the subleading logarithms coefficients are under control in the resummation formal- 
ism. At the same time, the NLL predictions for partonic cross sections begin to grow rapidly with the order of the 
perturbative expansion in due to the renormalon ambiguities |^ 10|. For this reason, the perturbative expansion 
is usually stopped at NNLO, avoiding the theoretical problem with power corrections. Nevertheless, we shall see 
in the next section that in the case of single spin asymmetry it makes sense to analyze the higher order corrections 
too. 



IV. HIGHER ORDER CONTRIBUTIONS 



Formally, the finite-order expansion procedure can be extended at NLL level to arbitrary higher order in as 



According to (2.14) and (2.15|), the LL contribution can be written as 



d^7s, , T"^ 3/2 , j^. , ^ /«,C^A"2n 
-^{s,^)=j^ dyxl''^—{xyS,ip)U{y) + ^ 



1 2n-l T 

In y 



(4.1) 



where y = s^/m?, Xy — (l — ym? / s^ , y-m = {s — 2m^fs)lvr? and 



da. 
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diy9 



2^ 



-2/3 + 2 In 



1+13 
1-/3 



-2/3(l-/32) + (l-/34)lni±| 



(1 + 7't,cos2(^) ^ (4.2) 



with /3 = -^/l — 4m2/s. At hadronic level, we have 



(4.3) 
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where 



(A)a51 iS) = i 



a 

with (f)g/ii{z) the gluon distribution function, and the quantity 



In* 2/ 



(4.4) 



-I + 



MS) = 



(4.5) 



describes the partial contribution of the ith power of the threshold logarithm. 

We have calculated numerically the quantities Ai{S) up to i = 11, i.e., up to the 6th order (N^LO) of perturbative 
expansion in ag. The results of our computations are presented in Fig. 6. One can see that Ai{S) are close to the 
LO result, ABor„(5) = Aa^-" (5) 7^^°™ (S), for all i < 11. 
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FIG. 6. Energy behavior of the Ai{E^)-ia,ctoTS in b- and c-quark production. 



So, we see that, at higher orders of perturbation theory, the LL predictions for the azimuthal asymmetry are also 
practically insensitive to soft radiation. The same situation takes place at the NLL level too. Indeed, neglecting 
in (2.14)-( 2T^ ) the small contributions of order of (1 — 2Cf/Ca) = l/-^c ^-^d the terms proportional to ln(ii/ui) 
which vanish after integration over ti and ui, we can write 

A^^^^ {S) ^ ± (l + {n 1)^) A,^MSy!^r'' (S) + O {llNj) , (4.6) 



with the partial contributions Ai{S) defined by (4^). Moreover, one can conclude from Eqs.(0)-(4^) and Fig. 6 
that all those subleading logarithmic (SLL) contributions which have the form 



(A)a^;:^ iS) 



C{n,i) 



S ry,, 

ds(l)g/h{s/S) I 

4m2 Jo 



dy4/2(A)a^;" {xy, 



y 



(4.7) 



-I + 



with arbitrary numerical coefficients C{n,i), also can not affect essentially the LO predictions for the azimuthal 
asymmetry. 

We have also investigated how the soft-gluon corrections to the asymmetry depend on the choice of the gluon 
distribution function. Parametrizing the gluon density as z~^(l — z)^, we have variated the parameters A and L 
in wide intervals: 0<A<1, 0<L<6. The energy behavior of the functions Ai{S) was found to be different 
at extreme values of L and A however, at L > 2, the ratio Ai{S) / A-BomiS) is equal to unity with a good accuracy 
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irrespective of A. Moreover, the softer the gluon density (i.e., the larger L) the smaller the corrections to the 
asymmetry. All curves given in Fig. 6 are calculated at A = 1, L = 5.|^ 



To clarify the origin of perturbative stability of the asymmetry, let us rewrite the Eqs.(4.1 )-(4.4) in a more 



convenient way. Changing the order of integrations in (4.4), we obtain 

UsCa \ 1 



(A)a^L(^) 



E 

n=0 



dz(^g/„(z)(A)<7^-" (zS) ^n{z,S), 



where 



^3/2 



At 
I372 



(f>g/h{T) In^" [{t - ^/zt) S/n? 



(4.8) 



(4.9) 



— lim 



In^ 



In^ 



eS_ 

2m? 



In^"-! [(7 



i7) 5/m2 



Formally, we can write 



E 



2n 



OisCA 



e5 



_20g/^(z) (z/4 + e + ^{z/4 + 0z/^ 



- 1 



1 

-<^>n{{A)Zn{S),S) dz^g/h 



(4.10) 



An? / S and grow at large z ~ 1 



with the mean value points Am? j S < z„(S'), Azn{S) < 1. 

At energies up to S/Am? ^ 10, the functions <i>„(z, 5), n> 1, vary slowly at z 
as In^" (1 — y/z). Since the gluon density rapidly vanishes at large z, (j)g/ii{z) ^ (1 — z)^, both values Zn^S) and 
Az„(S') are of the order of Am? / S. In other words, sufficiently soft gluon distribution function makes the collinear 
gluon radiation effectively soft at hadron level. In this case, the different high energy behavior of the spin-dependent 
and unpolarized Born cross sections is irrelevant since practically whole contribution to the radiative corrections 
to these quantities originates from the threshold region. Furthermore, a^™^{s) and Aa^™^{s) take their maximal 
values practically at the same energies: one can see from Figs.l and 2 that both Born level cross sections have 
peaks at 7y ~ 1. According to the saddle point arguments, these properties of the gluon distribution function and 
photon-gluon fusion provide with a good accuracy the equality^ 



$„ (A2„(5),^)~$„(z„(5),5), 



(4.11) 



which leads to spin- independent radiative factor in (4.10). 

As noted in previous Section and shown in Fig. 3, soft-gluon corrections to the partonic cross section depend 
essentially on the azimuthal angle (p. So, the mere spin-independent structure of the Sudakov logarithms can not 
explain our results. Our analysis shows that two more factors are responsible for high perturbative stability of 
the hadron level asymmetry. First, at the energies under consideration, the gluon distribution function supports 
the contribution of the threshold region. Second, the extremum point of the Born level cross section, 77 ~ 1, is 
practically 93- independent. 



■^The gluon distribution function is presently unknown beyond the NLO. We have verified that, in the case of CTEQ5M 
parametrization, Ai{S ) / Aso-n-iiS) is equal to unity to within 4% for all i < 11. 

^Note that Eq.(4.11) does not hold at very high energies since the functions 'bn{z,S) vary strongly in whole interval 
Am^/S < 2 < 1 at S/Arn^ > 10'\ In this case, the different high energy behavior of the spin-dependent and unpolarized Born 
cross sections becomes important and leads to (^-dependent soft-gluon corrections. 
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V. CONCLUSION 



In this paper we have investigated the impact of the soft gluon radiation on the single spin asymmetry in heavy flavor 
production by linearly polarized photons. Our calculations of the spin-dependent cross section at the NLL level up 
to the 6th order in as show that the azimuthal asymmetry is practically insensitive to the soft-gluon corrections at 
fixed target energies. Taking into account the remarkable properties of A{S) observed in our previous paper p^ , we 
conclude that, unlike the unpolarized cross sections, the single spin asymmetry in heavy flavor photoproduction is an 
observable quantitatively well defined in pQCD: it is stable both parametricaly and perturbatively, and insensitive 
to nonperturbative contributions. This asymmetry is of leading twist and can be measured at SLAG where the 
linearly polarized photon beam will be available soon [ p^Jl5| . Measurements of the azimuthal asymmetry would 
provide an ideal test of the pQCD applicability to heavy fiavor production. 
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